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Abstract. We consider the signed density of the extremal points of (two- 
dimensional) scalar fields with a Gaussian distribution. We assign a positive 
unit charge to the maxima and minima of the function and a negative one to its 
saddles. At first, we compute the average density for a field in half-space with 
Dirichlet boundary conditions. Then we calculate the charge-charge correlation 
function (without boundary). We apply the general results to random waves and 
random surfaces. Furthermore, we find a generating functional for the two-point 
function. Its Legendre transform is the integral over the scalar curvature of a 
four-dimensional Riemannian manifold. 
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1. Introduction 

In two-dimensional systems with orientational degrees of freedom, it is often fruitful to 
focus on the nodal points (zeros, defects), and to ignore other degrees of freedom. The 
distribution of the nodal points keeps information about the degree of (orientational) 
order within the system, and about the topology of the 'substrate'. Well known 
examples are defects in liquid crystal films, vortices in thin superconductors and 
superfluid films (for an overview, see pQ), or chaotic wave functions in quantum 
mechanics 0, microwave billiards [HI El, and optical speckle pattern |HJ. 

Here we take a closer look at the distribution of the nodal points of a gradient 
field or likewise the extremal points (maxima, minima, saddles) of a scalar field <j) with 
a Gaussian distribution^. The locations of the extrema of 4> are the points, where <fi 
is stationary, i.e. the gradient V(f> vanishes. The distribution of the stationary points 
was discussed earlier in the context of random surfaces [5] [7| |S] and of excursion sets 



The extremal points carry a topological charge q, which is the sign of the Jacobian 
right at the extremum q = sign(det(<9i<9j</>)) = ±1. In two dimensions, the positive 
extrema are the maxima and minima of the function </>, the negative extrema are the 
saddle points. We will compute moments of the signed density 



where r a are the positions of the extrema with sign (charge) q a . 

X The distribution of </> is solely characterized by the correlation function (0(r)0'(r')). Higher 
moments are calculated with the help of the Wick theorem. 
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From the experimental point of view, measuring the absolute mean values of the 
density, \p\, is much easier than measuring the usually small density difference between 
the density of the positive points minus the density of the negative ones. The signed 
density, however, is the much more appealing object from the mathematical point of 
view. The signed density is subject to a topological constraint, and its averages are 
differential geometrical objects, as shown in jllj . 

We represent the signed density of the extremal points of a d-dimensional scalar 
field <j> through [HE] 

p(r) = det {didj<P(r)) S d (V0(r)) . (2) 

In order to calculate the average of the absolute density (\p\), one needs the 
simultaneous distribution of the d(d + l)/2 random variables didj(f> and of the d 
variables V</>, where both sets are Gaussian random variables. Therefore, we need 
to know all moments (didj(f>dkdi(f>), (didj(f>dk<j>) , and {di<f>dj<j)) at coinciding points. 
For the averaged signed density, however, we need less information. According to 
it is sufficient to know the moments 

9ij {r) = {d^drfir)) (3) 

(which usually depend on the position r). The mean density (p) turns out to be 
proportional to the (total) curvature of the d-dimensional Riemannian manifold which 
is described by the above metric tensor g^. In two dimensions, the mean density 
is proportional to the Gaussian curvature K of that manifold ('surface') (times the 
invariant area element ^/det gtj ) 

2n(p(r)) =K^tg~. (4) 

Also, higher-order correlation functions fit into this scheme. The correlation function 
(p(ri) . . . p(r/)) is the total curvature of a / x d dimensional manifold with metric 
tensor § 

gia,jp = {di^(r a )dj(l)(r^)) , i,j = l...d,a,(3 = l...f (5) 

where a pair of a Greek and a Roman index is seen as a single composite index. The 
charge-charge correlation function for a Gaussian field in two dimensions is therefore 
the curvature of a particular four-dimensional Riemann manifold. We shall see, that 
the 'Einstein'-action, i.e. the covariant integral over the scalar curvature R of the four- 
dimensional manifold, plays an important role as a generating functional, relating the 
two-point function and the correlation function of the scalar field <fi. Before we discuss 
the two-point function, we present a simple example of a field <p with a nonzero mean 
density (p(r)), namely a system with a straight boundary. The plots of the charge- 
charge density and of the charge density near a wall are somewhat similar, since the 
latter quantity might be seen as the correlation of a charge with its mirror charge. 
The general results are applied to random waves, and to thermally fluctuating surfaces 
(fluid membranes). 

§ We define a / X d dimensional Gaussian scalar function <J>(ri, . . . , rj) = 4 } ( r i) + • ■ • + 4>( r f)- It i s 
obviously p$(r±, . . . , rj) = p(r±) X ... X p(rj), where p<j is the density of extrema of the field <!>. We 
obtain the metric tensor (diaQdjpQ) = (di(f>(r a )dj4>{i'p)) . 
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2. Density of extremal points near a boundary 

We calculate the signed density of extremal points of a Gaussian distributed random 
function </>(x, y) in the upper half space y > and assume Dirichlet boundary 
conditions <fi(x,y = 0) = (see Qlj). We begin with a field (f>(r) in Euclidean, two- 
dimensional space with an isotropic and translationally invariant correlation function 

y)<f>(x',y')) = G (y/(x - x' f + (y - y') 2 ) • (6) 

The correlation function for the upper half space is obtained by subtracting the 
appropriate mirror image 

G+(x, y\x, y') = (<j»(x, y)(f>(x', y')) 

= G ^/{x-x'f + {y-y'f) - G (y/(x - x'f + (y + y') 2 ) . (7) 

The Green function obeys the boundary condition G + (x,y\x' ,y' = 0) = 0. The 
components of the metric tensor (|3J) are easily obtained 

9ij(x,y) = (d i cf)(x,y)d j <f>{x,y)) : 

9xx = - G (0) H , 

2y 

g yy = - G"(0) - G"(2y), 

9xy = 9yx = 0. (8) 

The signed density of extremal points is the scalar curvature of a two-dimensional 
manifold with metric tensor JH1 



Anp(y) = y/fcU)R (9) 

where the scalar curvature is twice the Gaussian curvature R = 2K in two dimensions. 
The metric is orthogonal and depends on the distance to the wall y only. The non- 
vanishing affine connections are 

(d x d y (/)d x (j>) = d y g xx /2 
(dydy^d y 4>) = d y g yy /2 

(d x d x (/)dy4>) = - d y g xx /2 (10) 

since (d x (f>dy(p) = and, therefore, {d x d x cj)d y <j)) J r{d y d x (f)d x (j)) = 0. The scalar curvature 
is according to JT] 

R = (det g^eijBki^didk^djdicp) - (did h <t>d m 4>) g mn {d n (j)djdi4>)) 

7;d y g xx + \d y g xx {{gyy)' 1 d v g vv + {g^^dyg^)) (11) 



gxxgyy \ 2 ^ 4 
where e±2 = 1, e2i = — 1, en = e22 = 0. A straightforward calculation yields 

4Trp(y) = -dy {[gxxgyy)' 1 ' 2 d y g xx ^j = d y f{y) (12) 
where f(y) = - (g 

xxgyy ) d y g xx is the integrated charge density. 
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2.1. Random waves in half space 

Chaotic wave functions are well modelled by a superposition of partial waves with 
Gaussian distributed amplitudes and fixed energy El- The partial waves (j>k must obey 
the Helmholtz equation — V 2 0fc = k 2 <f>k with k being fixed. Without lost of generality 
we can set k 2 = 1. For time reversal symmetry, i.e. in the absence of a magnetic field, 
we can choose real amplitudes. The corresponding free-space correlation function of 
the cj) reads [TJ El HI ESI HI Ej 

G ( r ) = T^vT fd 2 P 6(p 2 - l)exp(ipr) cx J (r), (13) 



(27T) 2 , 

where Jo is the zeroth-order Bessel function. G(r) also models the correlation of the 
undulations 4>{r) of capillary waves (with fixed energy) at the surface of a liquid. 
We might therefore also apply the following results to the extremal points of random 
capillary waves or even to ocean waves El El E] ■ 

We compute the density of extrema of a random wave near a straight wall with 
Dirichlet boundary conditions. We expect it to be a model for the distribution of 
extremal points of a chaotic wave function near the boundary |12| . According to 
section El we have to compute the nonvanishing components of the metric tensor (PJ 

g xx = 1/2 - '^pl = (1 - J (2y) - J 2 (2y)) /2, 

g yy = 1/2 + J[(2y) = (1 + J (2y) - J 2 (2y)) /2 (14) 

which can be found in El- We obtain d y g xx — 2J2(2y)/(2y) and the integrated charge 
density (El 

M = =^ = — =^ -Tf 2 = -d + .78 +...)• (15) 

V9**9yy 2 y((l- J 2 (2y)) 2 -(J (2y)) 2 ) 

Figuredshows f(y) and its derivative 2np = d y f(y)/2. Next to the wall one observes 
a negative density, i.e. an excess of saddle points, followed by a sharp positive peak. 
The system has, contrary to the free one, a net excess charge 4ir J dyp = /(0) = 1/2 
per unit length. 



2. 2. A surface of revolution 

There is a good way to visualize the charge density i|12|) . It is possible to embed 
the manifold with metric tensor JSJl in three-dimensional space. This allows us to 
literally see the defect density, which is the curvature of the embedded surface times the 
element of area Vdet g. There is, however, a price to be paid - we must compactify the 
dimension x by making it 27r-periodic. Then one can construct a surface of revolution 
with metric tensor (JjJJ . The surface of revolution is parametrized through a coordinate 
y along the principal axis and the angle < x < 2n 

I A(y)cosx \ 

X(x,y)= A(y) sin x (16) 
V B(y) J 

with A > and a monotonically growing B(y). X is a vector in three-dimensional, 
Euclidean space. For a compact overview of the differential geometry of surfaces, 
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Figure 1. Integrated charge density /({/) (solid) and charge density 2irp (dashed) 
of a chaotic wave function next to a wall. 



especially of embedded surfaces, see e.g. 
surface of revolution reads 



j^El]. The (induced) metric tensor of the 



g xx =d x X-d x X = A 2 {y), 

9xy = 0, 

Qyy = 8yX ■ d y X = {dyAf + (dy Elf . 

Thus, the above metric is equal to the metric J3J) provided we identify 
My) = V<h^, 

d y B(y) = (g yy - (d y Aff 2 = Vfe^fp 



(17) 



B(y) = / dyy/g yy (y)(l-p(y)/4). 



(18) 

It is important to keep the information about the coordinate y, e.g. by displaying 
a grid with lines at equidistant intervals in y-space to be able to reconstruct the 
element of area y/det g. The embedding Ijl6|l of the metric tensor (PJ) is restricted to 
regions, where the integrated charge density obeys |/(y)| < 2. The contour of the 
surface of revolution for the case of random waves is shown in figure [21 where the 
function A(B) for < y < 7 is displayed. The crosses display the y-coordinate- 
two consecutive crosses represent an interval Ay = 1/8. Figures 03 and 01 show 
two projections of the three-dimensional surface. The meridians are at values of 
y = fc/4, k = 0, 1, 2, . . . , < y < 7. To obtain 4np for a given y, one has to find the 
meridian of number 4 x y (counted from the cone) , the local Gaussian curvature there 
and multiply it with ^/detg, which is proportional to the area of the corresponding 
grid element. The surface of revolution has a conical form for y — > 0, which is an 
artefact of the embedding (|16(l . In fact, expanding the surface for small y we obtain 
a conical geometry: A(y) — y/2 + . . . and B(y) = \/3y/2 + . . .. The curvature itself 
is negative at the origin and has the limiting value R = f'(y)(gxxg V y)^ 1 ^ 2 — * —1/2. 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 

Figure 2. Contour of the surface of revolution. Shows A(B) for < y < 7. Two 
consecutive crosses correspond to an interval of Ay = 1/8. 
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Figure 3. Surface of revolution. Shown for < y < 7. Two consecutive meridian 
correspond to an interval of Ay =1/4. 
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Figure 4. Different view on the surface of revolution. Shown for < y < 7. Two 
consecutive meridian correspond to an interval of Ay = 1/4. 




Figure 5. Scalar curvature R(y) for random waves. 



Figure displays the curvature as a function of the coordinate R = R(y). Note the 
pronounced positive peak of the scalar curvature at y s» 2. 



The distribution of extremal points of Gaussian scalar fields 



8 



2.3. Pinned fluctuating surfaces 

As a further application, we calculate the density of extrema for a nearly planar, 
fluctuating fluid membrane, attached to a straight line y — (imagine a membrane 
suspended in a large frame). We represent the shape of the membrane by a height 
variable <f){r). The domain of r is restricted to the upper half space y > 0, the 
boundary condition is (f>(x, y = 0) = 0. The thermal shape fluctuations of the fluid 
membrane are described by the Helfrich-Hamiltonian |2."{l 1241 

H/T = j d 2 r (k(V 2 0) 2 + cr(V</>) 2 ) (19) 

plus higher-order terms, k is the so-called bending rigidity, a is the effective surface 
tension. For almost planar surfaces the fluctuations of the height variable <j> are 
Gaussian to lowest order. We rescale the height y/nc/) — > <f> and set r = ajn. The 
correlation function for the unconstrained 4> is now 

G< "> = Wf ! d2p 7$TT) =h { ° + bs(1/r) - ' (20) 

where C is an irrelevant constant. We rescale the length scale that r = 1 
and regularize G at distances r < a. With the help of K' (r) = —K\{f) and 
K i( r ) = -(!/ 2 )(^o(r) +K 2 {r)) and K\{r)/r = -(l/2)(K Q (r) - K 2 (r)) we obtain 
the components of the metric tensor (jHJ 

9xx — B ~ (2y)- 2 + K 2 (2y)/2 - K (2y)/2 

g yy =B- (2y)- 2 + K 2 (2y)/2 + K (2y)/2, (21) 

where B = — G"(0). For distances y much larger than the correlation length (which 
is 1 in our units) we have g xx = g yy w B — (2y)~ 2 . The corresponding charge density 
<|12fl has an algebraic decay 
3 

4ttp - — y~ A for r > 1, (22) 
2B 

a reminiscence of the long ranged interactions (from the l/p 2 -term in the Green 
function). The integrated charge density f(y) is 

= -l/(2y 3 )+K 2 (2y)/y = -l/(2y 3 ) + K (2y)/y + K^yj/y 2 ^ 

a/ 9xx9 y y \/9xx9 y y 

For small arguments a < ?/ <C 1 the integrated charge density behaves as f(y) ~ 
y _1 F(log?/) + 0(y~ 2 ). Figure shows yf(y). Figure displays the function 4ny 2 p, 
where again y -2 is the leading singular behaviour of the charge density for small 
arguments. The cutoff region < y < a is not shown in both figures. 



3. The two-point function 

We present a convenient way of calculating the two-point function (charge-charge 
correlation function) C(r) = (p(O)p(r)). We have mentioned in the introduction 
that the two-point function can be expressed through the Ricmannian curvature 
tensor Rkr/,i\,ia,ji3 °f a four-dimensional Riemann manifold with metric tensor 
(di(f>(r a )dj<f)(rp)), where a = 1,2, f3 = 1,2 enumerate the two points. The Riemann 
curvature tensor is evaluated in |Appcndix B| The derivation of the correlation function 
with the help of differential geometry, however, is rather involved. Instead we utilize 
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Figure 6. The function yf(y) for membranes. 




Figure 7. The function 4-iry 2 p for membranes. 



a less tedious approach, which is even valid for more general Gaussian vector fields 
Ui,i = 1,2 with an isotropic, translationally invariant correlation function 

(u i (r)u j (r , ))=Xi j (r-r'). (24) 
In | Appendix A| we obtain for the charge-charge correlation function 

(2vr) 2 (p(r)p(O)) = (2vr) 2 C(r) = -$8^ (25) 
where the potential reads 



Qim{r) = eije mn -(deth) 1/2 e a be c ddjXacd n Xbd 



(26) 
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with hij — (v,iU m ) (ujU m ) —XimXjm- Summation over double indices is implied. Next, 
we restrict the vector field Uj to gradient fields Uj = &i4>. The nodal points of «j, i.e. 
the points where Ui vanishes, are the extremal point of the scalar field <p. The potential 
becomes 

toim(r) = eije mn ^(Aethy 1/2 e ab e cd djd a d c Gd n d b d d G (27) 

and hij = <5y — did m Gd m djG. For convenience, we have chosen the amplitude of <f> 
that — G"(0) = 1. Therefore, it is —didjG(Q) — 6ij. As a consequence of the polar 
invariance, matrices (such as ddG, f2 and h) have the general form 

'•< r < , a.. f x .. m. 



A a = A " + ^2 [Sit - -^-) . (28) 

For the particular choice r = (r, 0) the only remaining components are An and ^22- 
The non-vanishing derivatives of order 2, 3, 4 of G(r) and components of h and £ for 
r = (r, 0) are given in equation (|29|l (the prime ' denotes the derivative with respect 
to r) 

d x diG =Z y = G"(r) 
d 2 d 2 G =Z 2 = G\r)/r 
dxdxdtG = Z[ 
d 1 8 2 d 2 G =Z' 2 = (Zi - Z 2 )/r 
dxdxdxdxG = Z'{ 

dxdxd 2 d 2 G = Z'i (29) 
d 2 8 2 d 2 d 2 G = 3Z 2 /r 
hxx =Dx = l- Z\ 

h 22 =D 2 = l- Z\ 

£11 = Zx/Dx 

£22 = Z 2 /D 2 

A short calculation yields 

Q 1X = -(ZtfiDxD^- 1 / 2 , il 22 = Z' 1 Z' 2 (DxD 2 )- 1 /* (30) 

where Z x = G"(r), Z 2 = G'/r, D x = \- Z\,D 2 = \-Z\. We find for a matrix of the 
form (J2HJ 

didjAij = W 2 Axx + X -d r (Axx - A 22 ) , (31) 

where V 2 Ax\ = (rA' xl )' Jr. The correlation function becomes 

(27r) 2 C(r) = -didjSkj = -V 2 tt n - (fi n - fl 22 )' Jr = tf/r, (32) 
with the 'potential' 

1>(r) = n 22 - fi n - rQ'xx = "22 - (rO u )' 

3(Z' 1 -Z! l )(Z 1 -Z 2 ) (Zx - Z 2 ) 2 / 1/2 n ' 

- r( Dl D 2 )^ + r \ (DlD2) 

1 ({Zx-Z 2 f 



r(Zx - Z 2 ) V [DxD 2 y/i 



(33) 
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This result was also obtained by |26| using a direct approach. The correlation function 
of the Gaussian field <f> can now be written as 

G(r) = J d Pl dp 2 G (^JpI+pIJ cos(pir) = J dpG(p) cos(pr) (34) 

with G{p) > and also G{p) > 0. It is 

1 = -G"(r = 0) = / ApG{p)p 2 , (35) 



i.e. w(p) = p 2 G(jp) has the properties of a probability density. The r-expansion of 
G(r) reads 

„. , r 2 6r 4 cr e dr s er 10 „. 10l , „ 

G ^ = ~T + H 6f + ~8! ioT + °( r )' < 36 > 

where b = G""(r = 0) = J dpw(p)p 2 and c = — G^(r = 0) = J dpw(p)p 4 (as well as 
d and e) are positive constants. With the help of Mathematica™ we find the following 
series expansion of ip(r): 

4 6 (b 2 - c) r 2 



^-373' 3V3 

(456 4 -566 2 c + 3c 2 + 106d) r 4 
1 V ; +0(r 6 ). (37) 



540 V3 b 

It is shown in equation (|40|l that — ip(0) is related to the absolute density of extrema 
uq = (\p\). The next term of the tp expansion determines the charge-charge correlation 
function at almost coinciding points|| 

(2tt) 2 lim C{r) = ~ 2(C ~ &2) 

dpw(p) P 4 - ( / dpw(p)p 2 ) I . (38) 



3V3 \. 

The latter expression is obviously negative-the vicinity of an extremum is populated 
by extrema of the opposite type [23I2E|> which screen the central charge. 

3.1. Perfect screening of gradient fields 

We show that the charge-charge correlation function C(r) obeys the first Stillinger- 
Lovett sum rule 

lim / d 2 rC*(r) = -n . (39) 
£ ^ a J\r\>e 

It expresses the fact that a particular extremum is completely screened by extrema 

(charges) of the opposite sign jTHl 123 ED E2 Ej ■ Applying the results of 

section |3 we find for the left hand side 

roc 1 r°° 1 2(7"" (V = 0) 

-n = 2n rdrC(r) = — dr/(r) = - — #)) = r— 1 (40) 

Jo+ 27r Jo+ 37TV3 

j| C(0) is not well defined, since C(r) has a <5-likc singularity at the origin, see [HI . The (5-function, 
however, is lost 'somewhere' during the course of sectional A careful analysis at the origin similar 
to the one presented in 1111 is needed to recover the 5-function. 



The distribution of extremal points of Gaussian scalar fields 



12 



where G(r) is the correlation function of the field 4>. On the other hand, the absolute 
density of extremal points of the Gaussian field (f> reads [7| |Sj (for an arbitrary point 
r) 

no = {\d x d x cl>d y d y <t>-{d x d y <t>f\5{d x 4>)5{d y (t>)) . (41) 

The Gaussian variables di<j) are independent of the second derivatives didjtfi (evaluated 
at the same point). It is therefore no = (\ . . . |) (5(d x (j))5(dy(j))} with 



(6(d x <P)S(d y( f>)) 



d 2 E 



1 



1 

1/(270, 



/ E 2 



V(W) 2 



S(E x )6(E y ) 



^<(V0)2) (42) 

where we have chosen — G"(0) = 1. The second moments of didj<f> at coinciding points 
are 

{d x d v cj)d x d x 4>) = {d x dy4>dyd v 4>) = 

(d X d X (f>d X d X Ct>) = {dydydpdydy(j)) = G""(0) 

{d x d y <f>d x d y cl>) = (d X d x( f>dydy<f>) = G""(0)/3. (43) 

Next, we define the mutually independent Gaussian variables X = (d x d x (f> + 
d v dy^)/V2, Y = (d x d x <j>-d y d y 4>)/V2, and Z = d x d y <t> with (X 2 ) = (4/3)G""(0), 
(Y 2 ) = (2/3)G""(0), and (Z 2 ) = (1/3)G""(0). Then 

(|...|) = (\X 2 /2-Y 2 /2- Z 2 \) 

' ^ dxdydz exp (-(a; 2 + y 2 + z 2 )/2) 



(2tt) 3 / 2 



2g ,w (0) 2 G^(0) 2 2 

— 5 — x o — (y + z ) 



(44) 



We introduce spherical coordinates y = r sin# cos <p, z = rsin0sin<p, and x = rcosO, 
yielding 

rTr/2 



= G""(0) 



3V27T 



dr r exp(— r /2) 



d6> sin 6> 1 3 (cos (9) 2 - 1| 



4G W, (0) 
3\/3 



Finally, we obtain the absolute density of extremal points of < 
2G""(0) 



n 



3ttV3 

in agreement with equation l|4U|l . 



(45) 



(46) 



3. 2. Generating Junctionals for the two-point function 

We have found a representation of the two-point function C{r) as a functional 
derivative. Although this section is less useful from the practical point of view, it 
highlights an important mathematical aspect of the theory and demonstrates once 
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more the applicability of the differential geometrical approach. We begin with the 
integral over the two-point function times the correlation function of the field 0, 



-- b.c. + (27T)- 1 J A 2 rQ. ij d i d j G 
J rdr{Q 11 Z 1 + n 22 z 2 ) 
= b.c. + j rdr * {Z[Z' 2 Z 2 - Z' 2 Z' 2 Z X ) 

J dr {[D 2 /D 1 ) 1 / 2 Z[ + (D 1 /D 2 ) 1 / 2 Z^ 



= b 



= b 



c. 



n- - I dr-^—^iZi + Z'), 



= H-£ (47) 
where b.c. represents contributions from the boundaries r — 0,r — > oo, and 

D 

(Di 

C = - Jdr ({D 2 /D l ) 1 / 2 Z' l + (Di/Da) 1 /^^) (48) 

with D = 1 — Z\Z 2 . The action TC is identified in | Appendix B| as the integral over the 
scalar curvature^ R = g k "i< %a g lx iiP Rk lt i\ t i a ,jp 

H = {Ait)- 1 J d 2 rVdet gR. (49) 

This functional is well known from general relativity. The curved spacetime manifolds 
which obey the least-action principle 57i = are the solutions of the Einstein field 
equations (without matter). We show below that action l|49|l is also meaningful for 
our manifold. 

The action C is the generating functional of C(r) as computed in |Appendix C| 

\^=^/r={2-KfC{r). (50) 
The Legendre transform of C reads, using l|47() and the above equation 



The 'Einstein'-action H, expressed as a functional of the charge-charge correlation 
function C(r) is therefore the generating functional of the conjugated field-field 
correlation function G(r), 

1 SH[C] 



(27r) 2 r 5C(r) 



-G(r). (52) 



K In principle, one has to integrate over four coordinates. The metric, however, depends on the 
distance \r A — r B \ only. The integration over the centre of mass is therefore trivial. 
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3. 3. Two point function for random waves 

We now compute as an application the two-point function C (r) for Gaussian random 
waves <f> (without boundary). We rescale the correlation function (|13fl of the field 
4> for practical reasons and use G(r) = (</>(r)(/>(0)) = 2Jo(r). It is —G"(0) = 1, 
Zi = G"{r) = -Mr) + J 2 (r), Z 2 = -{Mr) + Mr)), Di = 1 - (J (r) - J 2 (r)) 2 , 
D 2 = 1 - {M r ) + M r )) 2 , yielding for the potential (J3HJ) 

*M = "ittt" ( uf- m) ■ ( 53 ) 

rJ 2 (r)dr \(1 - (J - J 2 )^ 2 (1 - (J + J 2 )^ 2 J 

With the help of the relations Jq = — J\, = J\ — 2J 2 /r and 2J\/r = Jo + J 2 we 
obtain the result 

2(J 2 ) 3 



ip(r) 



(1 - (Jo - J 2 ) 2 ) 1/2 (1 - (Jo + J 2 ) 2 ) 1/2 
1+ M , 2(J 2 - J ) 



•h) \.h ' 1-(J 2 - Jo) 2 , 
£ /_3_ J 2 - Jo J 2 + Jo 

r 2 U l-(Jr Jo) 2 1 - (J 2 + Jo) 2 



(54) 



The small argument expansion of tp reads 



1 / r 2 r 4 IWr" 



ib(r) = - 7 =\l )+0(r b ). (55) 

V3\ 48 2304 29859840 J 

The small r behaviour of the charge-charge correlation function (|32|l is therefore 

(2n) 2 C(r) = l — — + +0(r 6 ). (56) 

24\/3 V 24 207360 J y ' 1 ; 

Figure ISI shows tp(r)/2 and (27r) 2 C(r) for < r < 10. The plateau of the correlation 
function for small radii is almost perfect due to the smallness of the higher order 
expansion coefficients. Note the pronounced, negative peak at r » 3.4. 



4. Conclusion 



We have studied the signed density of extremal points of scalar field with a Gaussian 
distribution. We have calculated the average signed density in half space with Dirichlet 
boundary conditions. The density is expressed as the Gaussian curvature of an 
associated (abstract) manifold, depending on the correlations of the scalar field. The 
abstract manifold could be visualized by embedding it into three-dimensional space. 
Then we computed the two-point correlation function of the above density, which 
is in principle the total curvature of a four-dimensional manifold. We identified the 
'Einstein'-action as the generating functional, which relates the correlation function of 
the scalar field and the correlation function of the density of extrema. Currently, no 
application of this equation is known. Nevertheless, we hope that its further study will 
ultimately lead to a deeper understanding of the differential geometry of the abstract 
manifold, which governs the distribution of the extremal points. 
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Figure 8. The charge-charge correlation function (2ir) 2 C(r) (dashed) and its 
potential tf)(r) (solid; if)/2 is shown for scaling reasons). 



Appendix A. Computation of the general charge-charge correlation 
function 

We calculate the charge-charge correlation function for Gaussian vector fields u.i,i = 
1,2 with an isotropic, translationally invariant correlation function 

(«<(r)u J -(r')>=X«(r-r'). (A.l) 
The charge density p is represented as 

1 a ( eije k iu k djUi\ 1 Q _ a e 



2tt \ u 2 + e J 2 ir(u 2 + e) 2 

-> det^UjO^W for e ^ (A.2) 

where ei2 = 1, e2i = — l,en = e22 = 0. The charge-charge correlation function 
C(r) = (p(r)p(r')} is 

4A W = /e !jefcl e„e s ; tW ^ [(r) ^ /)9 « H ' (/) 



3r 4 \ " J ** st ( M 2 (r) + e)(w 2 (r') + e) 
= - did m Q im (r - r') (A. 3) 

utilizing the translations! invariance of the correlation function (|A.1|> . With the help 



of the Fourier-transformation 

2vr 



(u 2 + e) 1 = — I d 2 pexp(ipu)K (pV^) (A.4) 



The distribution of extremal points of Gaussian scalar fields 



16 



{K v is the second modified Bessel function of order v) we obtain 

Sl im {r) = J d 2 pd 2 qK (y/ep)K (y/eq)e lj e k ie mn e s 



d d d d 



dp k dq s dAji dB n 



a=b=o 



x (exp (AijdiUj(Q) + BijdiUj(r) + ipu(0) + iqu(r))) 

1 f, 2 ,2 PkQs d d 

a pa qeije k ie mn e s t- 



A=B=0 



x (exp (AijdiUj(0) + B t jdiUj(r) + ipuifS) + iqu(r))) (A.5) 

in the limit e — > 0. The latter expression was obtained with the help of two partial 
integrations and diKo(y/ep) = —\feK\{y/ep)pi/p — > —pi/p 2 . The Gaussian average 
can now be done, yielding 



ihm{r) = - — 7 / d pa qeije k ie mn e s t 



A=B=0 



x exp (-{piPj + qiQj) (uiUj) /2 - piqj (ui(0)uj(r))) 

x exp(AijBki (diUj(0)dkUt(r)) +ip l B a (u*(0)9jtuj(r))) 

x exp (iA i: jq k (diUj(0)u k (r)) + Terms cx A 2 ,B 2 ) 

1 f ,2 A 2 Pk Qs 

~4^J P 1^ kl e mn e st ^ 

x ((djUi(0)d n u t (r)) - p v (u v (0)d n u t {r)) q w (djUi(0)u w (r))) 
x exp (-(PiPj + qiQj) (uiUj) /2 - p*qj (uj(O)uj(r)}) 
1 



4^2 / d pd qe tj e k ie rnn e si ^ 2 ^ 2 

x {-djd n xit{r) + p v dnXvt(r)q w djXiw(r)) 

x exp (-(p^ + (tiitij) /2 - p^Xij (r)) , (A.6) 

where (uiUj) — (ui(0)uj(0)) = (m(r)uj(r)). We eliminate the second derivate djd n xit 
by a partial 'integration' 

O im (r) = e mn d n (^-^ J d 2 pd 2 qe ij e k ie st ^^d : jXit(r)exp(. . (A. 7) 

+ 4tt2 ^ ^P^^ ei i ekiemnest ^T ( ^ q T' ( d jXltdnXvw ~ dnXvtdjXlw) 

x exp (-(piPj + qiqf) (uiUj) /2 - Piq^Xtj) (A. 8) 

The term (|A.7|) does not yield a contribution to the charge-charge correlation function 
oc -did m Vt im since it is a pure curl. The term (d 3 xitd n Xvw ~ djXlwdnXvt) is 
antisymmetric with respect to the indices t, w: 

djXitd„Xvw - djXlwdnX bdjXladnXvb- (A. 9) 

The q 2 -denominator is cancelled, since et w e s tqsqw = —q 2 , yielding 

1 f j2 j2 



^im(r) = - ^ J a pd qe lj e kl e mn —p-e ab d : jX.iadnXvb 

x exp [-(piPj + q t q 3 ) (uiUj) /2 - PiqjXij) , (A. 10) 

where we have omitted the curl (|A.7|) . Now f2j m is no longer symmetric with 
respect to i,m; nevertheless only the symmetric part of fli m contributes to the 
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charge-charge correlation function. Equivalently, we might replace e a bdjXi a dnXvb by 
the symmetrized version e. a b{pjXla,®nXvb + dnXiadjXvb)/2- The latter expression is 
antisymmetric in l,v and can therefore be written as ei v e a be c ddjXca,dnXdb/2- Up to 
an antisymmetric component, Qi m reads 

fl im (r) = eije mn -^ J d 2 pd 2 qe ab e c ddjXcadnXdb 

x exp {-(piPj + qiq-j) (uiUj) /2 - piqjXij) ( A -U) 
Next we perform the p, q-integral, which is Gaussian now, and obtain finally 

^im(r) = eije mn ^(det hy 1/2 e ab e cd djXacdnXbd, ( A -12) 

where hij = (uiU m ) (ujU m ) — XimXjm- As an examination, we use equations (|A.3|I and 
(|A.12() to rederive the well known two-point function for the special case of vector fields 
with independent components Xij ( r ) = o~ijG(r) in two dimensions (see [TTIIT^IISTII^) ). 
Then 

O iTO = e.ije mn ^ djGd n G = eije mn djK(r)d n K(r) (A. 13) 

where K(r) = arcsin(G(r)/G(0)). Plugged into (|A.1|) we get 

47r 2 C(r) = -e^emndjdraKd^K = 2det(d i d j K) = 2K"(r)K'(r)/r (A.14) 
in agreement with the known result. 

Appendix B. The 'Einstein'-action 

We will now calculate the scalar curvature R and the 'Einstein'-action H = 
(47r) _1 J d 2 ry / det gR for a four-dimensional Riemannian manifold with metric giajp — 
—didjG(\r a — rP\). In principle, one has to integrate over 4 coordinates. The metric, 
however, depends on the distance \r A — r B \ only. The integration over the centre of 
mass is therefore trivial. We represent the correlation function of the field <j> through 
its Fourier transform 

((/)(r A )(t>(r B )) = G(\r A -r B \) = f d 2 pG{p) exp (ip ■ (r A - r B )) (B.l) 



The 4x4 metric tensor reads (£, j 6 {1, 2}; a, (3 £ {A, B}) 

dip* i 

>Q r ia Q r j r 

where 



ftoj/, = Mr 8 ^)) = / d^lt, (B.2) 



i,(r A ,r B ;p)^^G(p) £ exp(zp • r a ) (B.3) 

a=A,B 

is an immersion of the four-dimensional manifold in function space. Explicitly, we 
have 

9ia,j/3 = SapSij - c af3 didjG{r) (B.4) 
where c a p = 1 — 5 a p. The inverse metric tensor is easily calculated with the help of 

g ia ' j0 = (S a ph ij + c a/s &) (B.5) 
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where fty = dij — did m Gd m djG, — h^ 1 ^ and £ ik = didjGh? k . The determinant 
of the metric tensor is detg = deth. The components of the affine connection are 

r fc7;4QJ 73 = J d 2 pd ia d j ptp*d k j'ip = 5 a pe a ~,didjdkG(r), (B.6) 

where the index a is not summed. The components of the Riemannian curvature 
tensor are calculated from 



Rkn,ix,ia,jp = j d 2 p{d ia dk^*djfjdi\i}j) 



d z P {d ia d ki rd m ^) g m ^ nv J dV {d nv rd j0 dix*l>) 
- (ia) <-> (jf3). (B.7) 

Explicitly we have 

J d 2 p (diadk^djfjdaip) = 8 ai 8 l 3xc a pd l d :j d k diG{r) + 8 a f il sZ tlk i (B.8) 
and 

Rk~/,lX,ia,jP — (dafdfjX ~ SaxSff^Ca^didjdkdiG 

- S aj e a ^did k d m G g m ^ nv 5 xepvdjdid n G 

+ & fh e fill .d j d k d m Gg m i i ' r »'5 a xe aa ,d i d l d n G (B.9) 

where a and (3 are not summed, Zijki — didjd k diG\ r =o and (5 Q( g 7 5 = lifa = /3 = 7 = (5 
and zero otherwise. The Riemannian curvature tensor can be decomposed into 

Rk~/,ix,ia,ji3 = e a pe 1 x Sij k i - Sap-fxeijCkiu) + c a pCjxeijeki@/2, (B.10) 

where we have used (5 a7 5px+5p 7 5 a x)c a p = c a pCy\ and {S^Spx-Sp^ax^ap = e a pe 7 x 
(no summation over a, f3) and the definitions 

Sijki = didjdkdiG + (didkdpGdjdidqG + djd k d p Gdidid q G) £ pq /2 
to = e l3 e kl d l d k d m Gd :j did n Gh mn /2 

9 = eijekididkdmGdjdidnG^n/2. (B.ll) 

We obtain easily 

g k ~t' ia R k7 lX ia jl3 = (6pxh lk - cpx^k) Sijki + 5i3xeijekih lk (-uj + 9/2) 

+ cpxeijekitikQ/2 (B.12) 

and 

R = g k ^ ia g a ^R kj ,ix, ia ,j0 

= 2 - Mi) Sm - ^ + ^ + edete) . (B.13) 
To simplify this expression, we define 

Sijki = Tijki - {28i 3 5ki - Si k Sji - 6u6jk)Q/2, (B-14) 
where T^ki can be written as 

T^ki = didjdkdiG + didjd p Gd k did q G^ pq . (B.15) 
The scalar curvature is now 

R=2 (( h ** l -^ T «»-^ + £k)- (B.i6) 
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Furthermore 

h ik h^d i d j d m Gd k did n GUn + 29/ dct h = h^did^Gh^dkdldr.GUn- (B.17) 

For convenience we choose r A — r B = (r, 0). The non- vanishing partial derivatives 
of G(r) for that particular choice of r are listed in equation (|29(l . We have uj = 
Z[Z' 2 /Dx - {Z! 2 ) 2 /D 2 and 

u _ 1 (Z&Zx-Zz) Z' 2 {Z 1 -Z 2 )\_ ( D 



Vteth VDITT 2 V Dx D 2 J \^DTD~ 2 ) < ' B ' 18 ^ 

where D = 1 - ZiZ 2 . The integral / d 2 r(det /i)" 1 / 2 cj is therefore zero up to boundary 
contributions from r = and r — > oo. The integral over the scalar curvature 
H = (47r) _1 J d 2 r\/det gi? becomes (up to boundary contributions) 

U = (27T)- 1 / d 2 rVdcTI - Zik^t) T ijkl ) 



poo 

/ rdrCDiAO 1 / 2 
Jo 



Di L>i£>2 r-D 2 V-Di Do D 



Z'{ 2Z%D 3Z 2 ( Z[ Z' 2 \ 2 Zi 



2 



1 Di Z?i.D 2 rL> 2 L>? £> 2 -D : 



2 

l 2 #|T)- (B.19) 



Z?iL> 2 V 27 DiDl 
Z'{ and Z 2 ' are eliminated with the help of two partial integrations 



/■OO rjll 

/ rdr{ Dl D 2 fl^ 
Jo u i 



= „.., + , „,„„,„.,. M -*L + 3^-(^, . lU . a „ 



/■OO 7// n 

/•OO 

= b.c. + 2 / rdr(£>i-D 2 ) 1/2 
Jo 



^2 z x Z 2 Z 2 ( Z 2 ) 



2 



yielding 



W = b.c.+ / rdr(Z?!^2) 1/2 FT 



2Z 2 L> 3Z 2 



rD\D 2 rD 2 

3 Z^ Z2 3 ( ^2 ) ^ ^1 ^ ( ^2 ) ^ ^2 2 Z\ Z^ ^1 ^2 

D X D 2 + D X D 2 ThTjf + D\D 2 



b.c.+ [°° rdriD.D^/ 2 (-- §- 



2Z 2 £> 3Z 2 



rDiD2 rD 2 
Z'-^ Z2 Z2 3 ( Z2 ) ^ ^1 2 f ^2 ) ^ Z2 2 ^2 ^2 



D X D 2 D X D 2 D X D\ D\D 2 



+ ^27^ • (B.22) 
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We replace a single Z' 2 = [Z\ — Z 2 )/r in the last four terms of the integral 



3Z^(Z 1 -Z 2 )Z 1 2Z' 2 {Z 1 - Z 2 )Z 2 2Z' 1 (Z 1 - Z 2 )Z 2 



DiD 2 



DiD 2 2 



D\D 2 



B.23) 



Identity l|B.18|) allows us to simplify the last two terms of the above equation 

-f) ^fdrzJ*" 
\fTTjh, \Di D 2 ) J \VTJiTh 



= b.c. 



, 2Z',D 
dr- ~ 



We gather terms oc Z[ and Z 2 and gain the remarkably simple result 

D 



H = - dr 



(Zi + z 2 y 



\/DiD 2 

up to contributions from the boundary. 

Appendix C. Functional derivative of the Lagrangian 

We calculate the variation of 



(B.24) 



(B.25) 



C = - [dr (pa/DOVa^ + (D 1 /D 2 ) 1 ^ 2 Z 2 



(C.l) 



under an infinitesimal variation of the correlation function G(r) — > G(r) + n(r). We 
introduce the shorthand notations d\ = d/dZ\ and d 2 = d/dZ 2 and obtain 



6C 



dr 



(D 2 /D 1 ) 1 / 2 V '" + (D 1 /D 2 ) 1/2 (v'/r)' 



+ 
+ 



(d^/D^ 1 ' 2 ^' + d 2 (D 2 /D 1 ) 1 / 2 V '/r) Z[ 
(di(Dx/D 2 ) 1/2 ri" + d 2 (D 1 /D 2 ) 1/2 r]'/rj Z' 2 
b.c.+ J dr[{d l {D 2 /D l ) 1 ' 2 Z' x +d 2 (D 2 /D 1 ) 1 ' 2 Z' 2 
+ (di(Dx/D 2 ) 1/2 Z[ + d 2 (Di/D 2 ) 1/2 Z'^ rf/ 



(C.2) 



2 V 



dr 



d 1 (D 2 /D 1 ) 1 ^ 2 r 1 " + d 2 (D 2 /D 1 ) 1 / 2 r]'/r) Z[ 



r) Z'o 



(C.3) 



= b.c. 



dr 



(d 1 (D 1 /D 2 ) 1 / 2 - d 2 (D 2 /D 1 ) 1 / 2 ) (ZW - Z[ V '/r) 

(C.4) 



where 



Y = 



b.c. + / drY{Z' 2 if - Z'-tf/r) 

Z\ — z 2 

{B x D 2 fl 2 - 



(C.5) 
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The variation becomes 



5C = b.c. - J drn' ((Y Z' 2 )' + YZ'Jr) 



!>.(■.+ / dr^ {2Y{Z[ - Z' 2 ) + Y' (Zj. - Z 2 )) (C.6) 



or 

1 AC 19 / (y(z x - Z 2 ) 2 



r 5G r d r \ r{Z x - Z 2 ) 



= if} 1 jr = (2-K) 2 C(r). (C.7) 
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